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Introduction 


During  the  past  tw  decades,  space  physics  has  progressed  from  missions 
whose  goal  was  a  rudimentary  exploration  of  the  near-earth  magnetosphere  to 
the  present  stage  in  which  rather  detailed  modeling  and  understanding  of 
magnetospheric  plaasa  processes  has  emerged.  Nonetheless,  because  of  the  vast 
scale  distances  involved  within  the  magnetosphere,  it  has  been  a  very 
difficult  problem  to  probe  the  system,  concurrently,  at  enough  different 
points  to  truly  inderstand  the  complex  relationships  between  its  different 
parts . 

Understanding  just  the  'quiet'  or  '  equilibrium  state  of  the  magnetosphere 
has  been  a  challenge.  Even  more  difficult  has  been  the  problem  of  understand¬ 
ing  the  dynanic  behavior  of  the  magnetosphere.  Ihis  dynamic  aspect  of  the 
magnetosphere  may  be  effectively  discussed  in  terms  of  energy  input  from  the 
solar  wind  into  the  magnetosphere.  Such  excess  added  energy  causes  the  asg- 
netospheric  system  to  move  out  of  its  equilibrium  state  into  a  more  energetic 
state.  In  sane  cases  this  gives  rise  to  a  very  large  scale  disturbance  (the 
geomagnetic  storm)  idiieh  in  turn  causes  worldwide  effects.  Mich  more  fre¬ 
quently,  however,  disturbances  within  the  magnetosphere  tend  to  be  sanewhat 
more  localized  involving  the  regions  connecting  to  nightslde  auroral  field 
lines:  such  a  disturbance  is  termed  the  magnetospheric  substorm.  (See  the 
paper  by  McPherron  [1973.  and  papers  thereafter)  for  a  discussion  of  a 
phenomenological  model  of  substorms.) 

In  order  to  inderstand  better  the  nonequilibrium  behavior  of  the 
magnetosphere,  a  period  (July  2R-30,  1977)  was  chosen  for  intensive  study. 
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This  period  was  characterized  by  the  development  of  a  large  geomagnetic  storm 
and  also  by  the  occurrence  of  several  magnatospherlc  substorms  [Manka  et  al., 
1981].  In  addition  this  period  offered  the  advantage  that  there  were  a  total 
of  12  earth-orbiting  spacecraft  positioned  at  widely  separated  points 
immediately  upstream  and  throughout  the  magnetosphere  and  these  satellites 
provided  data  coverage  of  plasma  and  field  changes  associated  with  the 
geomagnetic  storm  and  sub  storms. 

In  order  to  exploit  fully  the  information  provided  by  such  a  wide  array  of 
spacecraft  probes,  an  effort  was  made  under  the  aegis  of  the  International 
Nagnetospherlc  Study  (IMS)  to  assemble  researchers  who  had  data  from  satel¬ 
lites  for  the  28-30  July  1977  time  period.  In  May  of  1979,  approximately  10 
scientists  with  Interest  in,  and  data  on,  the  high-energy  plasmas  of  the  mag¬ 
netosphere  met  at  the  National  Space  Science  Data  Center  (NSSDC)  located  at 
the  NASA  Goddard  Space  Flight  Center  in  Greenbelt,  Maryland.  In  a  workshop 
setting  called  CDAW  2.0  (Coordinated  Data  Analysis  Workshop  -  2),  the 
researchers  studied  data  Wiich  they  had  earlier  provided  to  a  central  computer 
facility.  This  central  computer  allowed  data  from  any  sensor  on  any  satellite 
to  be  directly  canpared  with  concurrently-eeqsired  data  from  any  other  sensor 
on  the  sane  or  any  other  satellite.  Thus,  tatlike  most  prior  space  research 
situations,  experimenters ,  modelers,  and  theorists  had  at  their  ianediate 
command  the  data  required  to  address  many  questions  about  magnetosph eric 
dynsnics.  Variations  of  plasma  conditions  throughout  much  of  the  magneto¬ 
sphere  could  thus  be  described  for  a  large  mmber  of  different  regions. 

There  are  two  distinct  facets  of,  or  reasons  for  studying,  energetic 
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particles  within  the  earth's  magnetosphere.  The  first  of  these  facets 
reflects  the  intrinsically  interesting  question  of  where,  and  how,  these 
particles  are  actually  produced,  say,  during  magnetospheric  substorms.  The 
second  facet  of  energetic  particle  studies  is  a  very  practical  and  pragaatlc 
one:  given  that  such  particles  exist  (l.e.,  that  they  can  be  observed)  how 
can  these  particles  be  used  as  tracers  or  probes  of  large-scale  magnet:  ric 

processes?  The  CDAW-2  subgroup  6  research  team  attempted  to  explore  ei  sf 
these  avenues  associated  with  energetic  particles. 

Tne  types  of  studies  carried  out  by  subgroup  6  were  basically  four 
number : 

(1)  Timing  and  morphology  of  particle  injections; 

(2)  Variation  of  particle  phase  space  densities; 

(3)  Measurement  of  boundary  motions  using  ion  (proton)  gradient 
anisotropies;  and 

(9)  Adiabatic  modeling  (with  increased  particle  flux  (i.e.,  injection), 
convection,  corotation,  and  drifts). 

In  the  following,  we  will  discuss  our  findings  derived  from  each  of  the 
above  lines  of  inquiry.  Our  Initial  research  efforts  were  concentrated  on  the 
1200  UT  substorm  of  29  July.  This  was  the  last  and  largest  (AE  *  1200  t)  of  a 
series  of  sub  storms  that  occurred  on  29  July  following  a  worldwide  SSC  that 
occurred  at  0027  UT  [King  et  al.,  1981;  Wilken  et  al , ,  19811.  We  concentrate 

hare  on  measurements  msde  at  geostationary  orbit  (6.6  Rg)  (4iere  a  total  of  six 
spacecraft  made  extensive  observations  of  the  energetic  particle  behavior. 
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Observations 


Figure  1  Is  a  geocentric  solar  eollptlo  projection  of  the  positions  of  the 
six  primary,  near-geostationary  satellites  used  In  the  present  study.  The 
ATS-6  and  1977-007  spacecraft  were  located  very  near  one  another  at  *  0300  LT. 
ATS -6  had  NOAA,  Aerospace,  and  TRW  energetic  particle,  UCLA  magnetometer,  and 
UNH  plasma  experiments  on  board,  while  77-007  had  Los  Alamos  energetic  parti¬ 
cle  sensors  on  board.  The  Los  Almnos-instrtmented  spacecraft  1976-059  st 
*0700  LT  was  bracketed  by  the  GOES-1  and  -2  satellites  which  carried  NOAA 
energetic  particle  and  magnetometer  instruments .  Finally,  the  European  Space 
Agency  satellite  GEOS-1  (1.3  <  r  <  8  R_)  carried  a  complete  complement  of 
plasma  and  field  measurement  instruments  and  was  located  near  apogee  at  *  1300 
LT. 

General  geomagnetic  activity  for  July  29-30,  1979  is  shown  in  Figure  2 
[see  also  Manka  et  al.,  1981],  The  upper  panel  shows  selected  high-latitude 
magnetometer  station  records,  while  the  second  panel  from  the  top  shows 
H-component  magnetograns  from  five  standard  auroral  zone  stations.  The  third 
panel  of  Figure  2  shows  mid-latitude  stations  from  several  geographic 
longitude  sectors.  The  bottom  panel  summarizes  auroral  electrojet  activity  in 
the  form  of  the  AE(5)  index,  i.e.,  the  index  derived  from  the  five  auroral 
zone  stations  shown  in  panel  2. 

Particularly  evident  in  Figure  2  are  the  storm  sudden  commencement  (SSC) , 
due  to  an  interplanetary  shock  wave  hitting  the  earth  at  0027  UT  on  29  July, 
[c.f.  King  et  al.  1981  and  Wllken  et  al.,  1981 1  and  the  rapid  storm  mainphase 
development  thereafter.  These  features  are  seen  clearly  in  the  midlatltude 
magnetograms  of  panel  3.  Also  quite  evident,  especially  in  the  plot  of  AE, 
are  the  generally  disturbed  auroral  zone  conditions  on  29  July  and  the  large 


GEOCENTRIC  SOLAR  ECLIPTIC  PROJECTION 


Fig.  1.  Positions  of  the  geostationary  and  near-geostationary 
(GEOS-1)  spacecraft  used  In  this  study.  The  nominal 
magnetopause  location  In  this  solar  ecliptic  projection 
is  also  shown. 
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Fig.  2.  Common  scale  magnetograms  and  AE  indices  for  July  28-30, 
1977. 
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substorm  (AL  >  1000  y)  st  *1200  UT.  As  previously  mentioned,  it  is  the  1200 
UT  suustorm  upon  tfiich  we  coneentrste  in  this  paper. 


Energetic  Psrtlcle  Behavior  st  0300  LT 

Fig  ire  3  shots  energetic  electron  Huxes  measured  by  instruments  on  the 
1977-007  spacecraft.  The  five  energy  channels  shown  are  a  representative 
sample  of  the  *  20  electron  channels  (E  .10  keV)  available  from  77-007  and/or 
ATS-6.  All  electron  channels  at  *  0300  LT  exhibit  roughly  the  same  sequence 
of  events  with  a  pronoineed  fl ux  decrease,  or  'dropout,'  commencing  at  *  1135 
UT  [see  Fennell  et  al..  1981].  The  nuxes  eventually  diminish  from  *1  to  3 
orders  of  magnitude  (depending  on  energy)  but,  as  is  especially  clear  in  the 
30  keV  channel,  the  measured  intensities  remain  nonzero  throughout  the 
dropout.  Hence,  it  is  concluded  that  the  geostationary  spacecraft  at  0300  LT 
entered  a  region  of  much  reduced  electron  flux,  but  they  did  not  mserge  into 
the  extremely  low  intensity  region  of  the  high  tail  lobes.  The  most  likely 
explanation  is,  therefore,  that  77-007  and  ATS  entered  the  high-latitude 
plasma  sheet  between  /into  and  1155.  In  the  northern  'horn'  of  the  plasma 
sheet  it  w->uld  be  expected  that  energetic  particle  fluxes  (prior  to  sUbstorm 
onset)  were  lower  than  in  the  outer  trapping  zone,  but  higher  than  in  the  tail 
lobes . 

After  the  flux  dropout,  the  electron  intensities  appeared  to  recover 
simultaneously  at  all  energy  levels  to  slightly  more  than  the  predropout 
values.  At  1200  ITT  there  was  a  large  increase  of  electron  flux  and  this 
Injection  corresponded  closely  to  the  sharp  negative  bay  onset  seen  at  College 
(cf.  figure  2).  Note  that  lack  of  energy  dispersion  between  the  several 
energy  channels  suggests  that  the  electron  'injection  front'  extended  as  far 
east  as  0300  LT. 

At  1205  UT,  another  substantial  flux  increase  or  injection  took  place. 

This  was  largest  and  most  evident  in  the  higher  (E  >  100  keV)  energy  channels. 
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Representative  differential  energetic  electron  flux 
profiles  for  the  period  1130-1300  0T  on  29  July  1977  at 
the  0300  LT  spacecraft  position. 


This  injection  spike  was  also  simultaneous  in  all  energy  ranges  (i.e.,  without 
energy  dispersion)  and  this  again  allows  the  conclusion  that  the  Injection 
region  extended  as  far  eastward  as  0300  LT. 

A  third  flux  injection  event  (with  some  evidence  of  energy  d-ispersion) 
occurred  at  </»  1208  UT.  Mote  further  that  after  1225  UT  the  drifting 
high-energy  electron  population  apparently  moved  azLmuthally  arovxid  the  earth 
and  once  again  passed  over  the  spacecraft . 

Given  this  observed  electron  behavior,  we  now  turn  to  the  energetic 
proton  flux  variations.  In  Figure  4,  several  representative  low-  and 
mid-energy  proton  (ion)  channels  (Tom  ATS-6  are  shown.  The  18-20  keV  channel 
is  from  the  U.  of  New  Hampshire  pi a an a  experiment  while  the  other  (bur 
channels  (33-150  keV)  are  (Tom  the  NOAA  energetic  particle  experiment. 

Prior  to  1200  UT,  the  energetic  protons  in  the  range  15-150  keV  exhibited 
behavior  very  similar  to  that  of  the  energetic  electrons  seen  in  Figure  3.  A 
pronounced  flux  dropout  was  seen  after  *  1135  UT,  but  at  least  for  particle 
energies  up  to  many  tens  of  keV  the  flux  dropout  was  not  total.  This  further 
suggests  passage  of  the  spacecraft  into  a  region  of  reduced,  but  nonzero, 
flux.  This  again  argues  that  the  spacecraft  entered  the  high-latitude  plasma 
sheet  where  the  presub  storm  particle  fluxes  were  intermediate  between  the 
outer  trapping  region  and  the  tail  lobe. 

Following  the  dropout,  (as  with  the  electrons)  a  recovery  characterized  by 
several  complex  (lux  variations  was  seen  in  the  protons.  Mote,  however,  that 
the  lowest  energy  proton  channels  showed  little  evidence  of  pronounced 
Injections  of  •new*  particles  since  the  average  intensity  level  was  the  sane 
both  before  and  after  the  dropout,  ty  contrast,  the  higher  energy  proton 
channels  (above  *  50  keV)  appeared  to  show  a  recovery  to  approximately 
predropout  flux  values  (1155-1200  UT)  and  then  showed  large  flux  increases  at 
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Representative  differential  energetic  proton  flux 
profiles  for  the  period  1130-1300  UT  on  29  July  1977  at 
the  0300  LT  spacecraft  position. 


substorm  onset,  i.e.  at  1200  UT.  This  oenpares  wall  wit'.,  the  electron 
injection  morphology  of  Figure  3. 


The  high  energy  proton  data  acquired  concurrently  at  77-007  are  shown  in 
Figure  5.  The  upper  panel  shows  that  145-340  keV  protons  also  exhibited  flux 
-dropout- recovery  sequence  just  as  the  electrons  and  low-energy  protons  did. 

The  recovery  fluxes  (1155-1200  UT)  were  a  factor  of  *  4  higher  than  prior  to 
the  dropout  at  11 30  UT.  In  the  lower  panel  of  Figure  5  it  is  seen  that  the 
very  high  energy  proton  (E  >  0.4  MeV)  fluxes  were  quite  low  before  1135  UT  (J 

O  i  i  e 

<  1  am-s“  -sr”  -keV"  )  and  were  near  background  during  the  dropout  period 
(1135-1155  UT). 

In  Figure  5  it  is  seen  that  the  injection  features  described  above  for 
electrons  and  mid-energy  protons  were  perhaps  even  more  evident  in  the 
high-energy  protons  measured  at  77-007.  Thus,  the  injection  spikes  commencing 
at  1200  UT  and  +  1205  UT  becmee  progressively  sharper  and  more  dlatinct  up  to 
at  least  several  hundred  keV. 

One  of  the  most  striking  aspects  of  the  data  in  the  lower  panel  of  Figure 
5  is  the  appearance  of  the  very  regular,  periodic  proton  drift-echo  pulses 
[cf.  Belian  et  al..  19?8  and  Baker  et  el. .  1979).  As  has  been  we 11- 
doe  imen  ted  in  the  literature,  these  high-energy  proton  pulses  are  injected 
into  the  outer  radiation  zone  at  substorm  onset  and  maintain  their  discrete 
identity  sufficiently  long  to  drift  azlmuthally  around  the  earth  maty  times. 

In  this  case,  it  is  seen  in  the  0. 8-1.0  MeV  channel,  for  exmaple,  that  at 
least  four  'echo*  pulses  were  recorded.  In  a  more  detailed  analysis  section 
below  we  will  return  to  the  information  provided  by  the  drift  echo  data. 

Energetic  Particle  Behavior  at  0700  LT 

Figure  6  is  the  0700  LT  counter  part  to  Figure  3,  i.e.,  it  shows  represent¬ 
ative  energetic  electron  channel  measurements  for  the  1130-1300  UT  period  on 
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Fig.  5.  High-energy  proton  data  at  the  0300  LT  spacecraft 
position  for  1130-1309  OT  on  29  July  1977. 


18 


SUBSTORM 
I  ONSET 


I 


I 


19 


es  Figure  3  for  the  0700  LT  epececreft  position 


29  July.  Note  that  except  for  the  lowest  energy  channel,  there  was  a  gradual 
flux  decline  between  1 1 30  and  r  1200  UT.  There  was,  however,  no  evidence  for 
the  aajor  flux  dropout  seen  in  the  midnight  sector  (as  revealed  by  77-007/ATS 
observations) . 

After  v  1205  UT  there  were  substantial  flux  Increases  In  all  eleotron 
energy  channels.  These  increases  were  gradual  in  character  with  apparent 
energy  dispersion  effects  [A-rnoldy  and  Chan.  1969].  These  observations  are 
consistent  with  the  electrons  being  Injected  over  a  broad  front  near  local 
midnight  (actually  extending  as  far  east  as  *  0300  LT)  and  subseqoently 
drifting  eastward  to  the  76-059 /GOES  location. 

Energetic  proton  data  from  spacecraft  1976-059  (0700  LT)  corresponding  to 
those  shown  in  Figure  5  are  presented  in  Figure  7.  In  that  figure,  the  upper 
panel  summarizes  the  145-390  keV  flux  variations  while  the  lower  panel 
summarizes  the  variations  of  the  very  energetic  component  (E  >  0.4  MeV) . 

As  was  seen  in  Flgtre  6  for  the  electrons,  the  proton  fluxes  shown  in 
Figure  7  also  exhibited  a  gradual  flux  decline  prior  to  1200  UT,  but  showed  no 
major  flux  dropout.  The  behavior  of  the  proton  fluxes  at  0700  LT  after  * 

1200  UT  was  highly  energy  dependent.  Up  to  *  250  keV,  the  proton  fluxes 
appeared  to  recover  gradually  and  indistinctly  with  some  possible  energy 
dispersion.  By  contrast,  the  >0.4  MeV  proton  population  exhibited  a  very 
dear  onset  with  considerable  energy  dlaperslon.  As  la  evident  from  the  lower 
panel  of  Figure  7,  the  high-energy  proton  behavior  was  of  the  clear  drift-echo 
character.  Careful  comparison  of  the  details  of  shape  and  timing  of  the 
pulses  in  Figure  7  with  those  of  Figure  5  shows  two  things: 

(1)  The  pulse  shapes  at  03  and  07  LT  were  remarkably  similar  in  width  and 
amplitude  for  any  given  energy  channel;  and 

(2)  An  identifiable  drift-echo  pulse  in  any  given  channel  at  07  LT 
appeared  slightly  before  the  sate  pulse  sppesred  at  03  LT. 
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Same  as  Figure  5  for  the  0700  LT  spacecraft  position 


It  ia  concluded  that  eaaentially  all  of  the  proton  results  aeen  at  the 
76-059 /GOES  location  can  be  accounted  for  by  an  injection  of  protona  near 


midnight  with  a  subsequent  westward  drift  completely  around  the  earth  to  the 
0 f  LT  position.  The  complex  recovery  behavior  of  150-250  keV  proton  fluxes 
prior  to  *  1225  UT  appears  to  have  been  primarily  adiabatic  and  will  be 
discussed  further  below.  The  width  of  the  proton  injection  regions  around 
local  midnight  are  progressively  broader  at  lower  and  lower  energies. 

Energetic  Particle  Behavior  at  1W)  LT 
Figures  8  and  9  show  the  electron  and  proton  fluxes,  respectively, 
measured  at  the  location  of  GE0S-1 .  The  data  are  shown  in  the  form  of  stacked 
energy  spectra  in  each  instance.  In  the  case  of  both  particle  species,  the 
|  lowest  energy  channels  show  a  gradual  flux  decrease  between  <r  1130  and  1200 

UT  followed  by  a  gradual  recovery.  Only  in  the  higher  energy  channels  (E  >  80 
keV)  was  the  recovery  very  sharp  or  dramatic.  The  flux  recoveries  in  both 
particle  species  showed  very  clear  energy  dispersion  with  the  recovery 
occvrrlng  first  in  the  higher  energy  channels.  This  feature  is  consistent 
with  substorm  injection  of  energetic  particles  (broadly)  near  midnight  with 
subsequent  drift  of  the  particles  to  the  1300  LT  position. 

Phase  Space  Density  Variations 

In  the  foregoing  section,  we  have  discussed  pronounced  flux  increases  in 
terms  of  injections.  That  is,  we  have  presiaed  that  the  flux  enhancsnents 
actually  corresponded  to  new  or  'fresh'  particles  transported  to,  or 
accelerated  in  the  vicinity  of,  geostationary  orbit.  In  order  to  confirm  this 
supposition,  we  have  evaluated  the  particle  distribution  functions  at  constant 
first  adiabatic  invariant  [Schulz  and  Lanzerottl.  197k]. 
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The  phase  space  density,  or  distribution  function,  of  a  particle  popula- 
tion  can  be  defined  in  terms  of  adiabatic  invariants  and  time: 

f(u,J.t>  «  J/p2  (1) 

Here,  u  is  the  first  adiabatic  Invariant  (magnetic  mcment)  ,  J  is  the  second 
Invariant,  t  is  time,  J  is  the  directional  differential  particle  flux,  and  p 
is  the  relativistic  manentum.  Taking  J  to  be  the  equatorial  perpendicular 
flux  and  J  £  o,  the  phase  space  density  at  constant  first  invariant  is  given 
by 

f(u.t)  =  J/2m0wB  (2) 

where  mQ  is  the  particle  rest  mass  and  B  is  the  tot  i  equatorial  magnetic 
field  strength. 

As  is  evident  from  Eq.  (2),  the  advantage  of  studying  the  phase  density  at 
constant  u  is  that  adiabatic  (magnetic  field)  variations  are  removed.  Thus 
true  particle  density  increases  or  decreases  are  revealed  and  sources  or  sinks 
of  particles  can  be  identified.  In  particular  this  analysis  can  reveal 
vAiether  or  not  new  particles  were  injected  in  the  1200  UT  sub  storm  on  29  July. 
Figure  10  shows  exemplary  spectra  which  were  obtained  at  various  times  for 
this  event  period.  The  panels  on  the  left  show  spectra  for  the  01  LT 
spacecraft  grouping,  while  the  panels  on  the  right  show  similar  data  for  the 
07  LT  grouping.  The  upper  panel  in  either  case  shows  j  for  electrons,  vtille 
the  lower  panels  show  J  for  protons. 

As  is  evident  from  Figure  10,  the  data  are  distributed  relatively 
accurately  according  to  a  simple  exponential  spectrun.  This  is  particularly 
true  below  </•  300  keV.  The  dashed  line  accompanying  each  set  of  data  is  the 
least-squares  fit  to  the  observed  particle  distribution  where  the  fit  is  given 
by 


25 


PROTONS  {cm2- s-sr-kcW)*1  ELECTRONS  (cm2  -s-sr-keVf1 


Fig.  10. 


Differential  energy  spectra  for  electrons  (upper  panels) 
and  protons  (lower  panels)  at  the  UT's  indicated  near 
each  data  set.  Spectra  taken  at  ~0300  LT  are  shown  on 
the  lefthand  side  while  spectra  taken  at  ~0700  LT  are 
shown  on  the  righthand  side.  Dashed  lines  show  least- 
squares  fits  to  the  data  points. 


3  = 


(3) 


2  1 

and  has  units  of  particles  (cm  -s-sr-keV)  .  Except  in  the  highest  proton 
energy  ranges  after  ^1220  UT  (where  drift-echo  effects  are  dominant),  the 
spectral  fits  provide  an  excellent  analytical  representation  of  the  observed 
spectra.  Our  procedure  in  the  present  analysis,  therefore,  has  been  to  fit 
(for  each  1-min  flux  average)  the  observed  energy  spectrum  to  obtain  K(t)  and 
EQ(t).  Given  these  fits,  we  thus  have  j(E,t)  to  be  used  in  Eq.  (2). 

The  other  required  information  for  phase  density  calculation,  is  the  total 
magnetic  field  strength.  In  Figure  11  we  show  the  values  of  |fl|  for  the  0300 
and  0700  LT  spacecraft  locations.  The  largest  variability,  as  might  have  been 
expected,  was  seen  in  the  nighttime  sector  at  the  ATS-6  location.  Because  one 
component  (Y)  of  the  ATS-6  magnetometer  [McPherron  et  al..  19751  was 
inoperative  at  the  time  of  these  measurements,  the  inferred  field  line 
direction  from  electron  anisotropy  data  at  77-007  was  used  to  complmnent  ATS 
field  data.  In  a  CDAW-2  algorithm  procedure  the  two  measured  ATS-6  field 
components  (X  and  Z)  and  the  field  line  direction  from  77-007  were  sufficient 
to  provide  the  total  field  vector,  B,  at  v<n300  LT. 

Combining  the  magnetic  field  data  of  Figure  11  with  energy  spectral  data 
for  each  minute  between  </>ll30  and  1300  IJT  gave  us  the  desired  phase  space 
densities  at  constant  u.  The  ranges  of  u-values  selected  for  investigation 
were  chosen  as  follows.  The  miniraun  and  maximum  kinetic  energies  of  electron 
and  protons  measured  on  any  of  the  six  observing  spacecraft  were  first 
considered.  The  minimum  energy  measured  was  E^j  *  10  keV  whereas  the  maximum 
energy  channel  from  which  useful  data  were  obtained  was  E  *  1.0  MeV.  The 
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Fig.  11.  Total  magnetic  field  strength  variations  for  the  1200  UT 
substorm  at  the  0300  and  0700  LT  grouping  positions 
studied  in  this  paper. 
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measured  range  of  B,  similarly  considered, was  Bmln  *  100  y  and  B^^  *  25 0  Y • 

Thus,  the  u  range  was 


(*») 


and  with  some  spectral  extrapolations ,  a  useful  range  of  u- values  in  this  case 
was  *  1-1000  HeV/gauss. 

Figure  12  shows  examples  of  the  phase  space  densities  for  electrons  at  v  * 
1,  10,  and  100  MeV/G.  The  most  evident  features  in  the  upper  panel 
(77-007/ATS  grouping)  were  the  following: 

(1)  Even  with  removal  of  adiabatic  effects,  the  flux  dropout  persists; 

(2)  The  phase  space  densities  at  constant  v  were  identical  before  the 
dropout  (<r  1130  UT)  and  after  the  dropout  (*  1155  UT); 

(3)  True  phase  space  density  increases  were  observed  for  all  magnetic 
manents  (energies)  after  1200  UT. 

The  points  above,  therefore,  demonstrate  that  in  a  broad  sector  near  local 
midnight  there  was  a  large  scale  boundary  motion  which  took  the  observing 
spacecraft  into  a  low  density  region  (i.e.,  across  a  spatial  discontinuity). 
This  thinning-like  event  clearly  preceded  the  substorm  onset.  Prior  to  the 
substorm  onset  the  midnight- sector  spacecraft  also  returned  to  a  predropout 
density  configuration  for  several  minutes  (1155-1200  UT);  this,  therefore, 
clearly  was  not  an  Injection  event.  At  *  1300  UT  a  clear  injection  of  new  or 
'fresh'  particles  occurred  for  all  magnetic  moments. 

The  lower  panel  of  Figure  12  shows  the  electron  density  variations  at  07 
LT.  Comparison  of  these  results  with  Figure  6  shows  that  at  this  location 
essentially  all  flux  variations  before  *  1205  UT  were  adiabatic.  Viewing  the 
phase  space  densities  in  this  region  of  the  magnetosphere  shows  essentially 
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Fig.  12.  Electron  phase  space  density  variations  (computed  as 
described  In  the  text)  for  the  1200  UT  substorm 
period.  Densities  at  constant  first  Invariant  values 
(P,  as  labeled)  are  plotted. 


flat  profiles  prior  to  1205,  a  density  dip  at  *  1205,  and  energy-dispersed 
density  increases  after  *  1206  UT,  consistent  with  injection  and  drift  froai 
the  west . 


Proton  phase  space  density  variations  are  shown  in  Figure  1?  for  a  •  1, 
-50,  and  300  MeV/G.  Can  pared  to  the  electron  data  of  Figure  12,  much  more 
variability  was  seen  in  the  proton  density  profiles.  This  in  part  represents 
statistical  variations  in  the  measured  fluxes  which  translate  into  variations 
of  K(t)  and  EQ(t)  in  Eq.  (3).  Nonetheless,  the  following  features  seem  to  be 
established  by  the  data: 

(1)  At  03  LT,  there  appeared  to  be  a  phase  space  density  increase  for 
very  low  u-values  between  *  1135  and  1150  UT  while  at  higher  v-values 
a  clear  dropout  was  seen; 

(2)  Substantial  injections  of  new  particles  were  seen  at  03  LT  for  u  >  10 

iT 

MeV/G  but  little  clear  evidence  exists  for  injection  of  new  protons 
with  low  u-values; 

(3)  At  07  LT,  there  may  have  been  some  significant  dips  and  peaks  before 
*  1210  UT,  particularly  at  high  v-values,but  the  most  substantial 
effects  occurred  after  *  1210  UT  as  protons  azimuthally  drifted 
westward  from  midnight  to  the  76-059 /GOES  location. 

Gradient  Anisotropy  Tnformatior 

By  exanining  flux  and  phase  space  density  variations  (particularly  at  the 
03  LT  position),  it  is  established  that  new  partlcles(up  to  several  hundred 
MeV/G)  appeared  at  synchronous  orbit  between  <r  1200  and  1210  UT  on  29  July.  A 
remaining  question  about  these  particles  is  where  the  particles  came  from. 

The  best  available  tool  for  exmninlng  the  question  of  the  general  source 
region  for  the  injected  hot  plasms  and  energetic  particles  is  provided  by  ion 
gradient  measurements.  Because  of  their  large  gyrorsdll,  10-1000  keV  protons 
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13.  Proton  phase  space  density  variations  for  the  1200  UT 
substorm  period.  Densities  at  constant  first  Invariant 
values  (P,  as  labeled)  are  plotted. 
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can  provide  good  Information  about  density  gradients  that  exist  within  a 
region  of  strong  radial  intensity  variations  or  within  an  injected  cloud  of 
plana  and  energetic  particles  tFrltz  and  Williams.  1979:  Williams  et  el.. 
1979;  Palmer  et  al..  1976;  Walker  et  al..  1976]. 

.  The  spacecraft  77-007  and  76-059  are  particularly  well  suited  for 
exninlng  ion  gradient  anisotropies.  The  reason  for  this  is  that  these 
spacecraft  spin  about  an  axis  that  points  continually  toward  the  center  of 
earth  and  proton  fluxes  are  measured  by  sensors  with  view  directions  that 
point  radially  outward  perpendicular  to  the  spin  axis.  Thus,  a  rather 
complete  scan  of  ion  fluxes  Is  obtained  on  each  ten- second  spacecraft  rotation 
in  both  the  east-west  sense  and  in  the  norttwsouth  sense.  Given  the  fact  that 
100-200  keV  protons  have  typical  gyroradii  of  several  hundred  km  U  0. 1  R£)  at 
synchronous  orbit,  one  can  probe  regions  far  removed  from  the  spacecraft  by 
the  gradient  anisotropy  technique. 

The  gradient  paraneters  are  computed  as  follows: 

=  (E  -  W)/(E  ♦  W) 

vtoere  E  is  the  proton  flux  (E^  >  1**5  keV)  measured  in  the  sector  with  the 
detector  looking  eastward  and  W  is  the  proton  flux  measured  looking  westward. 
Similarly, 

ANS  =  (M  -  S)/<1»  ♦  S> 

»*iere  N  is  the  north-looking  measured  flux  and  S  is  the  south- looking  measured 
flux.  Given  the  direction  of  the  magnetic  field  in  the  vicinity  of  the 
satellites,  and  using  the  sense  of  gyration  of  protons,  Agy  >  0  generally 
implies  a  higher  density  (flux)  inside  the  S/C,  whereas  A^w  <  0  implies  a 
higher  density  outside  the  S/C.  For  a  stretched  (taillike)  m^netic  field 
orientation  (as  distinguished  from  a  completely  dipolar  field)  one  also 
obtains  some  secondary  information  from  A^.  Similarly,  the  primary 
information  from  ANg  concerns  higher  flux  above  the  S/C  ( A^  >  0)  or  below  the 
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S/C  (Agg  <  0).  The  implications  of  various  kinds  of  gradient  anisotropies  are 
summarized  in  Table  1. 

We  only  present  end  A^  for  the  03  LT  position  here  slnoe  this  mbs  the 
primary  reg  ion  into  which  the  direct  particle  injection  mss  observed.  In 
order  to  give  s  sense  of  the  magnetic  field  orientstion  st  thst  location. 
Figure  14  shows  the  msgnetlc  field  line  meridional  tilt,  #g.  The  solid  line, 
for  reference,  is  the  field  tilt  at  0700  LT,  while  the  dashed  line  is  the 
value  of  at  0300  LT.  Note  that  in  a  dipole  magnetic  field,  would  be  the 
magnetic  dip  angle  (Bg  =  tan”1 (2tanx ) ) .  For  the  76-059/GOES  spacecraft  this 
means  the  dipolar  value  would  be  *  25°  while  for  77-007/ATS  the  dipolar  value 
would  be  *  10*.  An  extreme  taillike  (nondipolar)  magnetic  configuration, with 
the  field  lines  lying  nearly  parallel  to  the  magnetic  equatorial  plane,  is 
seen  at  03  LT  dirlng  the  flux  dropout.  This  again  seems  to  reinforce  our 
interpretation  that  a  large-scale  boundary  motion  took  place  during  the 
dropout  period.  It  also  strongly  suggests  that  the  spacecraft  entered  the 
high-latitude  plasma  sheet  where  very  tallllke  field  would  be  expected.  We 
note  that  the  appearance  of  this  tail like  field  topology  is  a  common  precursor 
to  stibstorm  onset  [HcPherron,  1970,  Baker  et  si..  19781  and  apparently 
indicates  m  extreme  stressing  of  the  outer  magnetosphere  prior  to  the 
sub storm  energy  release  at  1200  UT. 

Figure  15  shows  the  (upper  panel)  and  Agg  (lower  panel)  values 
calculated  from  the  77-007  energetic  proton  data  (E  >  145  keV) .  Looking  at 

and  Agg  together,  the  following  sequenoe  of  events  is  seen.  Between  1155 
snd  *1200,  i.e.  during  the  recovery  from  the  flux  dropout,  was  strongly 
positive.  This  suggests  that  the  higher  particle  density  was  inside  the 
spacecraft.  A^g  during  this  sane  period  was,  for  the  most  part,  strongly 
negative,  suggesting  a  high  particle  flux  below  the  spacecraft.  Sinoe  Figure 
14  showed  the  field  to  be  very  tallllke  during  this  period,  our  contention  of 


34 


Table  1.  Proton  Gradient  Anisotropy  Information 


Parameter  Value  High  Density  Location  High  Density  Location 

(Primary  Information)  (Secondary  Information) 


*EV 

> 

0 

Inside  S/C 

( Below-East) 

aew 

< 

0 

Outside  S/C 

(Above-West) 

ahs 

> 

0 

Above  S/C 

(Inside-East) 

ans 

< 

0 

Below  S/C 

(Outside-West) 

29  JULY  1977 


Fig.  14.  Field  line  dipolar  tilt  angle,  A  ,  inferred  (aa 
described  in  text)  from  the  energetic  electron 
anisotropy  symmetry  axis.  A  very  taillike  field 
(®,  >  60°)  is  observed  at  ~0300  LT. 


Fig.  15.  Proton  (E  >  145  keV)  gradient  anisotropy  data  for  the 
period  1155-1225  UT  at  the  0300  LT  position.  The  east- 
west  anisotropy  is  shown  in  the  upper  panel,  while  the 
concurrently  measured  north-south  anisotropy  is  shown  in 
the  lower  panel. 
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a  boundary  motion  during  the  dropout  with  the  high  flux  region  mowing 
earthward  and  equatorward,  is  fully  borne  out.  As  the  fluxes  recover,  the 
spacecraft  were  enveloped  from  inside  and  from  below. 


At  1200  UT,  Agy  went  strongly  negative.  This  period  corresponded 

precisely  to  the  first  energetic  particle  and  hot  plasma  injection  into 

synchronous  orbit.  The  character  of  A^  showed  that  the  Injected  particles 

cane  from  outside  the  spacecraft  location.  For  this  sane  period  A  _  was 

NS 

strongly  positive, showing  that  the  particles  also  generally  arrived  from  above 
the  S/C.  The  conclusion  is  therefore  unambiguous  in  this  case,  vit.  the 
injected  particles  arrive  at  6.6  Rg.  from  the  outside  and  from  above.  This 
very  likely  means  that  these  particles  filled  the  high-latitude  plasna  sheet 
and  that  these  filled  field  lines  then  collapsed  inward  over  the  spacecraft. 

After  the  leading  edge  of  the  particle  injection  passed  over  the 
spacecraft,  went  strongly  positive  and  A^«.  was  weakly  negative  (1202-1205 
ITT).  This  indicates  that  the  highest  particle  density,  after  the  injection, 
iv-  <  generally  inside  6.6  R£. 

A  second  particle  injection  occurred  (cf.  Figures  3  and  5)  it  r  1205  UT. 

Figure  15  shows  again  that  these  particles  cane  from  outside  6.6  R£  since  A£W 

♦ 

was  strongly  negative.  Note  in  Figure  14  that  B  was  more  nearly  dipolar  by 
1205  UT.  Therefore,  in  this  case  A^,,  became  only  weakly  positive  during  the 
injection.  It  is  concluded  with  considerable  certainty  that  the  1205  UT 
injection  of  energe  ic  particles  and  hot  plasma,  as  was  also  true  for  the  1200 
UT  case,  cane  froa  outside  of  synchronous  orbit. 

The  apparent  1203  UT  injection  of  particles  (see  Figures  3  and  5)  seemed 
to  show  energy  dispersion  effects,  consistent  with  the  interpretation  that  the 
injection  front  did  not  directly  extend  as  far  eastward  as  0300  LT.  Indeed,  a 
substantial  gradient  anisotropy  signature  of  this  injection  is  not  seen  in 
Figure  15. 
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A  composite  plot  of  the  >  1H5  keV  proton  flux  and  the  computed  value  of 
a£W  is  shown  in  Figure  16.  The  recovery  sequence  between  1155  and  1200  UT, 
the  flux  injection  at  1200  UT,  and  the  flux  injection  beginning  at  1205  UT  are 
all  particularly  evident  in  that  figure.  Minor  (but  statistically 
significant)  changes  in  A^  between  1212  and  1225  UT  sre  also  evident  as 
subsequent  snail  pulses  of  protons  drift  past  the  spacecraft,  approaching  from 
the  east  (A^  >  0)  and  receding  to  the  west  (A^  <  0). 

Drift-Echo  Timing  Information 

Proton  drift-echo  events  such  as  shown  above  in  Figures  5  and  f  can  be 
used  to  infer  times  and  locations  of  the  'centroids'  of  particle  injections 
fBellan  et  al..  1978].  As  illustrated  by  the  detailed  10-s  flux  averages  in 
Figure  17,  the  pulses  of  drifting  protons  show  evidence  of  basically  a  triple 
structure  in  each  pulse.  These  more  detailed  (10-s)  flux  values  have  been 
used  to  carefully  determine  the  time  of  the  'pesk  1',  'peak  2'  and  'peak  3' 
relative  flux  maxima  for  the  0.H-0. 5,  0.5-0. 6,  0.6-0.R,  and  0. 8-1.0  MeV 
channels  at  the  0700  and  0300  LT  positions.  The  local  time  of  the  observed 
peak  pulses  (modulo  36O 0 )  was  then  considered  versus  the  universal  time  of 
each  peak  pulse.  Two  drift-echo  pulses  could  be  clearly  discerned  in  each  of 
the  two  lower  energy  ranges,  while  three  pulse  echoes  were  seen  at  the  two 
higher  energies. 

Table  2  summarizes  the  LT  and  UT  data  points  inferred  from  the  high- 
resolution  drift-echo  data.  It  should  be  noted  that  discerning  the  individual 
relative  'peaks'  was  uncertain  when  the  pulses  overlap,  fin  the  other  hand, 
sane  of  the  peak  times,  as  might  even  be  evident  from  the  1-min  averages  of 
Figires  5  and  7,  are  quite  distinct  and  obvious.  Other  of  the  peak  times  had 
to  be  judged  from  relatively  subtle  inflections  in  the  flux  profiles. 

Overall,  the  data  points  in  Table  2  have  associated  UT  uncertainties  of 
approx imatel y  ♦  1  rain. 


39 


Fig.  17.  High  time  resolution  plot  for  the  0.4-0. 5  MeV  proton 

channel  on  1976-059  at  ~0700  LT.  The  triple  peak  time 
structure  of  the  pulse  Is  demonstrated. 


Table  2.  Proton  Drift-Echo  Pulse  Time 


Peak 

Energy 

(MeV) 

S/C 

76-059  i 

CUT) 

S/C 

77-007 

CUT) 

#1 

0.4-0. 5 

1212.0 

1228.0 

— 

1215.0 

1231.0 

— 

0.5-0. 6 

1211.0 

1226.0 

— 

1214.0 

1228.0 

— 

0. 6-0.8 

1209.5 

1221.5 

1233.5 

1211.5 

1222.5 

1236.0 

0.8-1 .0 

1209.0 

1219.5 

1232.0 

1210.4 

1221.5 

1233.5 

12 

0.4-0. 5 

1217.2 

1233.5 

— 

1220.5 

1236.5 

— 

0. 4-0.6 

1216.3 

1231.0 

— 

1219.0 

1233.0 

— 

0 . 6-0. 8 

1214.5 

1225.5 

1237.0 

1216.0 

1227.0 

1238.5 

0 . 8—1  •  0 

1213.0 

1223.0 

1235.0* 

1214.5 

1224.5 

1235.0 

#3 

0.4-0. 5 

1220.5 

1237.1 

— 

1223.5 

1239.5 

— 

0.5-0. 6 

1219.5 

1235.0 

— 

1222.5 

1237.0 

— 

0. 6-0.8 

1217.0 

1228.5 

1239.5 

1218.5 

1229.9 

1241.0 

0 • 8—1  •  0 

1216.0 

1225.5 

1236.0 

1217.6 

1227.0 

1237.0 

‘point 

not  used 

in  least- squares 

fit. 
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Results  for  the  several  selected  76-059  and  77-007  energy  channels  are 
plotted  in  Figire  18.  In  each  panel  we  separately  plot  data  for  each  of  the 
peak  1  through  peak  3  pulses.  The  paraaeter  ♦  is  equivalent  to  LT  (in 
degrees)  except  that  it  runs  clockwise  from  midnight  (in  the  same  sense  as 
proton  drifts)  rather  than  counterclockwise.  For  each  energy  channel  a 
least-squares  fit  through  the  data  points  is  shown. 

As  seen  by  Figure  18,  it  is  possible  to  arrive  at  an  Internally  consistent 
interpretation  of  all  of  the  high-energy  proton  data,  at  both  0700  LT  and  0300 
LT.  This  interpretation  is  that  there  were  three  high-energy  proton 
Injections  centered  in  the  post-midnight  region  and  these  injections  each 
exhibited  several  echoes  that  were  individually  seen  at  both  the  0700  and  0300 
local  times.  The  universal  times  of  the  injections  Inferred  from  Figire  18 

are:  peak  1  events.  *  1200  Iff ;  peak  2  events,  *  1205  Iff;  and  peak  3  events,  * 
1208  UT. 


Adiabatic  Modeling  Results 

A  major  underlying  theme  of  our  analysis  has  been  that  substorm  energetic 
particles  are  injected  in  the  nightslde  magnetosphere  and  that  these  particles 
subsequently  are  trapped  and  drift  to  positions  removed  from  the  injection 
site,  rtjch  of  the  foregoing  analysis  has  been  carried  out  within  this 
franework  and,  generally,  supports  such  an  Interpretation.  However,  in  order 
to  model  the  injection  and  drift  more  quantitatively  the  time-dependent 
convection  model  of  Smith  et  al.  [1979]  was  used. 

This  model  follows  the  motions  of  charged  particles  under  the  Influence  of 

the  geanagnetic  and  electric  fields.  A  Volland-Stern  type  of  convection 

2 

electric  field  (E  *  — ?  ♦  and  •  =  AR  sin  4)  and  a  dipole  magnetic  field  are 


assumed.  Here  •  is  the  electric  potential,  ♦  is  a  local  time  paraneter 
measured  frcra  local  midnight,  and  R  is  geocentric  radial  distance.  As  shown 


•w 


Fig.  18.  Local  time  ('4>')  versus  UT  plots  for  high-energy  drift- 
echo  pulses  seen  at  S/C  1977-007  and  1976-059.  As 
discussed  in  the  text,  the  Intersections  of  the 
manifolds  of  lines  in  each  panel  give  an  idea  of  the 
local  time  and  universal  time  of  the  proton  injection. 
The  small  inset  polar  plot  in  the  central  panel 
illustrates  the  S/C  locations  and  the  centroids  of 
proton  injection. 
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by  Snlth  et  al.  [19791,  the  tine-variation  in  the  electric  field  may  be 
characterized  by  the  geomagnetic  index  Kp  which  is  then  introduced  via  the 
paraneter  A. 

Although  this  large-scale  convection  model  has  been  quite  successful  in 
predicting  the  behavior  of  low-energy  charged  particles  during  storms  [c.f. 
Smith  et  al..  19791,  a  goal  of  the  CDAW-2  effort  was  to  test  the  model  for 
higher  energy  particle  injections.  Figure  19  illustrates  several  of  the 
results  for  "high-energy"  trajectory  simulations.  In  each  ease,  protons  • « 
with  u  =  1.0  keV/t  (100  MeV/G)  and  pitch  angle  s  90°  were  injected  at  a 
boundary  of  10  R^.  For  y  =  lOOMeV/C,  the  kinetic  energy  of  the  protons  at 
L=6.6  would  be  about  100  keV.  In  Figure  19  one-hour  increments  of  the 
trajectories  between  1200  and  1600  UT  on  29  December  are  displayed.  In  each 
Instance  the  GEOS-1  orbit  is  shown  for  reference. 

Case  A  shows  the  nominal  model  calculations  for  assumed  proton  injections 
at  2300  LT  through  0300  LT.  As  is  evident  from  the  figure,  the  normal 
convection  model  described  in  the  preceding  paragraph  gives  rise  to  untrapped 
particle  drift  trajectories  which  typically  encounter  the  dayslde  boundary 
near  local  noon.  In  Figure  19B,  the  Volland-Stern  convention  field  was 
decreased  to  one- fourth  strength  in  the  radial  range  £>-10RE  and  1000-1900  LT. 
This  change  causes  the  particles  to  be  "pulled"  back  on  the  dayslde  and  the 
relatively  high-energy  protons  injected  at  0200-0300LT  are  thereby  trapped. 
(Note,  however,  that  the  boundary  between  the  decreased  field  and  the  normal 
model  field  in  case  B  is  non-physical) . 

Magnetic  field  observations  in  the  outer  magnetosphere  during  the  substorm 
period  under  investigation  indicated  a  gradient  (aB/B)  value  much  less  than 
the  normal,  nonstorm  value.  In  case  C  of  Figure  19,  AB/B  was  reduced  to 
one-half  its  normal  value  to  be  more  consistent  with  observations.  This 
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feature  increased  the  trapping  efficiency  somewhat,  but  moat  trajectories  from 
the  midnight  sector  still  remained  untrapped. 

Finally,  in  case  D  the  field  gradient  was  maintained  as  in  case  C  but  the 
eonvection  electric  field  was  increased  by  a  factor  of  2  in  the  radial  range 
7-lORg  and  between  2000  LT  and  0400  LT.  This  change  drove  particles  more 
deeply  into  magnetosphere  initially  and  thereby  increased  the  trapping. 

(Again,  the  boundary  between  the  Increased  and  normal  field  is  unphysical.) 

In  suaunary,  the  time- dependent  convection  model  can  produce  trapped  drift 
trajectories  for  the  higher  energy  proton  component  (>100  keV) .  The  changes 
to  the  normal  model  in  order  to  accomplish  a  large  trapping  ratio  (such  as 
changing  the  magnetic  field  gradient)  appear  quite  consistent  with  observation 
and,  thus,  seem  to  provide  reasonable  physical  improvements  to  the  model.  In 
most  cases,  it  is  seen  that  only  high-energy  protons  injected  near  0200-0300 
LT  are  durably  trapped.  It  is  Interesting  that  our  proton  drift-echo  analyses 
also  tend  to  show  injection  positions  near  0200  LT  for  the  observed  proton 
pulses  in  this  sub  storm  case  (c.f.  Figure  IB). 


III.  DISCUSSION  AND  SUMMARY 


In  this  paper  we  have  used  energetic  particle  and  aagnetic  field  data  from 
six  satellites  near  geostationary  orbit  to  study  an  intense  substorm  period  on 
July  29,  1977.  Using  these  several  spacecraft,  well-distributed  in  local 
time,  has  given  us  a  perspective  on  global  substorm  phenanenology  not 
previously  available.  Several  different  analysis  techniques  (of  which  some 
are  unique  to  energetic  particles)  were  applied  to  the  data  sets  and  a 
self-consistent  picture  of  the  event  period  has  emerged. 

The  following  list  summarizes  our  observational  results  for  the  1200  UT 
substorm  at  the  three  local  times  sampled: 

Observations  at  0300  LT 

-  Taillike  magnetic  field  topology  was  seen  prior  to  substorm  onset. 

-  Large-scale  boundary  motion  occurred  as  indicated  by  the  flux 
dropout . 

Dropout  boundary  motion  was  to  the  inside  and  below  observing 
spacecraft . 

Observing  satellites  remained  in  a  finite  flux  region  (high-latitude 
plasna  sheet) . 

In  recovery  from  the  dropout,  the  spacecraft  were  enveloped  from 
below  and  inside . 

Two  clear  particle  injections  occurred  (1200  and  1205  UT)  with 
injection  fronts  extending  as  far  east  as  0300  UT. 

Injected  particles  clearly  came  from  outside  and  above  S/C. 

High-energy  proton  drift-echoes  were  seen  (injected  at  +  0100-0200 
LT) . 
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Observations  at  0700  LT 


-  Weak  flux  decline  was  observed. 

-  Only  mildly  taillike  magnetic  field  stretching  w&s  seen. 

-  Energy-dispersed,  injected  electron  population  was  observed: 

1205-1220  UT. 

-  Initial  proton  injection  spikes  were  only  weakly  manifested. 

-  Proton  drift-echo  peaks  were  clearly  seen  (injected  *  0100-0200  LT) . 
Observations  at  1300  LT 

Energy-dispersed,  injected  protons  and  electrons  (E  >  50  keV)  were 

\P 

observed:  1205-1220  UT. 

Most  low-energy  (E  <  50  keV)  particle  effects  (1130-1300  UT)  were 

ad iabatic . 

Based  on  the  results  presented  here,  some  very  firm  conclusions  regarding 
substorm  phenomenology  can  be  stated  and  these  results  can  be  extrapolated 
slightly  to  speculate  on  the  missing  pieces.  First,  there  seems  to  be 
considerable  evidence  that  the  magnetosphere  went  through  a  period  of 
substantial  energy  storage  prior  to  the  sudden  energy  release  at  *1200  UT 
[ McPherron ,  1970,  Baker  et  al..  1978].  An  attractive  and  consistent 
interpretation  is  that  this  energy  storage  manifested  Itself  as  a  taillike 
change  of  the  magnetic  topology  at  6.6  before  the  substorm  which  in  turn 

caused  the  observed  flux  dropout.  The  developing  magnetic  stress  seemed  to 
relax  slightly  (1155-1200  UT)  and  than  at  1200  UT  it  was  suddenly  relieved  in 
the  midnight  sector  simultaneous  with  the  injection  of  the  first  pulse  of  hot 
plaana  and  energetic  particles . 

Our  results  also  show  that  the  injected  substorm  particles  came  from 
outside  (and  above)  the  spacecraft  at  *  0300  LT.  With  the  present 
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information  we  are  uiable  to  tell  from  how  far  outside  6.6  R^  the  particles 
originated.  Given  the  very  stretched  magnetic  field  topology  that  existed 
during  the  injection  process,  it  is  quite  possible  that  the  field  lines 
carrying  the  injected  particles  actually  extended  deep  into  the  plasma  sheet 
(i.e.  beyond  10  R^).  One  point  that  is  clear  is  that  there  was  only  a  very 
low  level  of  energetic  protons  with  E  >  0.3  MeV  in  the  outer  radiation  zone 
before  the  substorm  onset,  and  yet  a  large  flux  of  such  particles  clearly 
appeared  at  geostationary  orbit  at  substorm  onset.  Adiabatic  modeling  shows 
that  trapping  can  be  simulated  by  convection  of  high-energy  particles  from 
beyond  lOR^. 

Several  recent  papers  have  discussed  the  outer  zone  plasma  injection 
process  in  terms  of  convection  electric  fields  [cf.  Kaye  and  Klvelson,  1979 
and  references  therein].  These  papers  show  that  inward  convection  of  plasma 
sheet  particles  associated  with  large-scale  substorm  electric  fields  can  lead 
to  substantial  particle  acceleration  (as.  Indeed,  was  the  case  in  the  modeling 
represented  in  Figure  19).  In  this  regard,  however,  Klvelson  [1980]  has  shown 
for  the  1200  UT  event  discussed  in  this  paper  that  acceleration  of  particles 

up  to  «r  1  MeV  cannot  be  done  with  the  usual  solar-wind  imposed  convection 
electric  field. 

Klvelson  [1980]  has  argued  that  the  sub  storm  induction  electric  field  may 
play  an  important  role  in  the  energization  of  the  high-energy  particles  seen 
in  this  event.  Using 


/•  100  y/S  min 


(see  Fig.  11) 


Kivelson  estimates  (using  V  x  f 
2 

A*  r  l  AB/At  , 

where  A*  is  the  change  of  electric  potential  and  t  is  the  scale  of  the  region 
in  which  B  was  collapsing.  Assuning  A*  is  of  the  order  of  1  MV  and  AB/At  ./>  20 
T/min ,  gives  t  *  9  R^.  Such  a  scale  size  for  the  region  of  near-tail  collapse 
associated  with  the  substorm  seens  reasonable  and,  thus,  suggests  that 
induction  fields  could  accouit  for  the  observed  particles  as  geostationary 
orbit . 

Based  on  large  numbers  of  other  high-energy  proton  events  observed  at 
synchronous  orbit  and  in  the  plasna  sheet,  Baker  et  al,  [1979]  argued  in 
favor  of  the  importance  of  induction  electric  fields.  They  showed  from  the 
timing  and  duration  of  energetic  proton  events  that  particles  with  energies  of 
*  1  MeV  cannot  be  produced  by  a  small  inward  radial  convection,  say  from  8-10 
Re;  large  impulsive  acceleration  must  be  responsible  for  their  production 
[e.g.  Pellinen  and  Heikkila,  19781.  The  high-energy  proton  results  shot*)  for 
this  event  are,  therefore,  consistent  with  the  plasna  sheet  energization  model 
presented  by  Baker  et  al.  [19791. 

In  sunnary,  it  sesns  evident  that  the  multiple-spacecraft  observational 
approach  used  here  is  a  powerful  one.  Since  the  geostationary  satellites  that 
we  have  used  in  this  paper  have  acquired  literally,  years  of  concurrent  data, 
we  look  forward  to  many  future  Joint  studies  of  the  effects  of  geomagnetic 
storms  and  substorms  on  magnetospheric  energetic  particle  populations. 


.  .  ||>  th.t 
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LABORATORY  OPERATIONS 


The  Laboratory  Operation*  of  The  Aaroapaca  Corporation  la  conducting  expar- 
lnantal  and  thaoratical  invest lgat Iona  nacaaaary  for  tha  evaluation  and  applica¬ 
tion  of  scientific  advances  to  new  military  apace  systems.  Versatility  and 
flexibility  have  bean  developed  to  a  high  degree  by  the  laboratory  personnel  in 
dealing  with  the  many  problems  encountered  In  the  nation's  rapidly  developing 
space  systems.  Expertise  In  the  latest  scientific  developments  is  vital  to  the 
accomplishment  of  tasks  related  to  these  problems.  The  laboratories  that  con¬ 
tribute  to  this  research  are: 

Aerophyslcs  Laboratory:  Launch  vehicle  and  reentry  aerodynamics  and  heat 
transfer,  propulsion  chemistry  and  fluid  mechanics,  structural  mechanics,  flight 
dynamics;  hlgh-temperaturc  thermomechanics,  gas  kinetics  and  radiation;  research 
In  environmental  chemistry  and  contamination;  cv  and  pulsed  chemical  laser 
development  Including  chemlcsl  kinetics,  spectroscopy,  optical  resonators  and 
beam  pointing,  atmospheric  propagation,  laser  effects  snd  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions,  atmo- 
spherlc  optics,  light  scattering,  state-specific  chemical  reactions  and  radia¬ 
tion  transport  In  rockst  plumes,  applied  laser  spectroscopy,  laser  chemistry, 
battery  electrochemistry,  space  vacuim  and  radiation  effects  on  matsrlals,  lu¬ 
brication  and  surface  phenomena,  thermionic  amission,  photosensitive  materials 
and  detectors,  atomic  frequency  standards,  and  bloenvlronmsntal  research  and 
monitoring. 

Electronics  Research  Laboratory:  Microelectronics,  GaAs  low-noise  and 
power  devices,  semiconductor  lasers,  electromagnetic  and  optical  propagation 
phenomena,  quantum  electronics,  laser  conmunlcetlons,  Ildar,  and  electro-optics; 
communication  sciences,  applied  electronics,  semiconductor  crystal  and  device 
physics,  radiometric  Imaging;  millimeter-wave  and  microwave  technology. 

Information  Sciences  Research  Office:  Program  verification,  program  trans¬ 
lation^ performance-sensitive  system  design,  distributed  architectures  for 
spaceborne  computers,  fault-tolerant  computer  systems,  artificial  intelligence, 
and  microelectronics  applications. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metal  matrix 
composites,  polymers"  end  new  forms  of  carbon;  component  failure  analysis  and 
reliability;  fracture  mechanics  and  stress  corrosion;  evaluation  of  matsrlals  In 
space  environment;  materials  performance  In  space  transportation  systems;  anal¬ 
ysis  of  systems  vulnerability  and  survivability  in  enemy-induced  environments. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radiation 
tTom  the  atmosphere,  density and  composition  of  the  upper  atmosphere,  aurorae 
and  alrglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation  of 
Plasma  oaves  in  the  magnetosphere;  solar  physics,  infrared  astronomy;  the 
effects  of  nuclear  explosions,  magnetic  storms,  and  solar  activity  on  the 
earth's  atmosphere,  ionosphere,  snd  magnetosphere;  the  effects  of  optical, 
electromagnetic,  and  particulate  radiations  In  space  on  epace  systems. 


